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Abstract. We report the first results of an ongoing spec- 
troscopic survey of blue horizontal branch stars in the 
metal-rich ([Fe/H] ~ -0.5) globular clusters NGC 6388 
and NGC 6441. Based on data obtained with the ESO- 
New Technology Telescope (NTT), we provide gravities 
and temperatures for four stars in NGC 6388 and three 
stars in NGC 6441. These results are marginally inconsis- 
tent with the predictions of canonical evolutionary theory, 
but disagree strongly with all non-canonical scenarios that 
explain the sloped horizontal branches seen in the colour- 
magnitude diagrams of these clusters. 



Key words: Stars: early-type - Stars: fundamental pa- 
rameters - Stars: horizontal-branch - globular clusters: 
individual: NGC 6388 - globular clusters: individual: 
NGC 6441 



1. Introduction 



The metal-rich globular clusters play a fundamental role in 
determining the formation history of our Galaxy (e.g. Min- 
niti |1995| , |1996j; Ort olani et al. |1995[ Zinn |1996[ Barbuy et 



al. 1998j ; Rich 1998 ). Moreover, they provide a crucial tem- 



plate for interpreting the integrated spectra of elliptical 
galaxies and for understanding the evolution of old metal- 
rich stellar systems. The horizontal branch (HB) stars in 
the metal-rich clusters are particularly important, since 
they provide a standard candle for determining distances 
(and hence ages) and are believed to be the major con- 
tributors to the UV- upturn phenomenon (e.g. Caloi 1989 ; 
Greggio & Renzini 1990, 1999; Bressan et al. |1994|; Dor- 
man et al 



1995; Yi et al. 1998) 
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Recent Hubble Space Telescope (HST) observations 
have found that the metal-rich globular clusters NGC 6388 
(C1732-447) and NGC 6441 (C1746 -370) ([Fe/H] = 
—0.60 and —0.53, respectively; Harris |l996|) contain an 
unexpected population of hot HB stars and therefore ex- 
hibit the well-known second parameter effect (Rich et al. 
1997 ). Most surprisingly, the mean HB luminosity at the 
top of the blue HB tail is roughly 0.5 mag brighter in V 
than the red HB "clump," which itself is strongly sloped 
as well. Differential reddening cannot be the cause of these 



sloped HB's (Piotto et al. |1997| ; Sweigart & Catelan [1998 
Layden et al. |1999[ ) 



The second parameter effect has often been attributed 
to differences in age or mass loss on the red giant branch 
(RGB). However, canonical HB simulations show that in- 
creasing the assumed age or RGB mass loss moves an HB 
star blueward in the V, B — V plane but does not increase 
its luminosity. Thus some other second parameter (s) must 
be causing the slope d HB' s in NGC 6388 and NGC 6441 
(Sweigart & Catelan [1998[ hereafter SC98). 

Three non-canonical scenarios have been suggested to 
explain the sloped HB's and long blue HB tails in these 
metal-rich globular clusters (SC98): 

1. High cluster helium abundance scenario: Red 

HB stars evolve along blue loops during most of their 
HB lifetime. For larger than "standard" helium abun- 
dances Y, these loops become considerably longer, 
reaching higher effective temperatures and deviating 
more in luminosity from the zero-age HB (ZAHB). 
If the cluster stars form with sufficiently high Y (Z 
0.4), the HB will slope upward (Catelan & de Fre- 
itas Pacheco 1996), as observed in NGC 6388 and 



NGC 6441. However, this scenario also predicts a much 
larger value for the number ratio R of HB to RGB 
stars than the value recently obtained by Layden et 



al. (1999) for NGC 6441. Thus a high primordial he- 



lium abundance seems to be ruled out as the cause of 
the sloped HB's and will not be considered further in 
this paper. 
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Rotation scenario: Rotation during the RGB phase 
can delay the helium flash, thereby increasing both 
the final helium-core mass and the amount of mass 
loss near the tip of the RGB. The net effect is to 
shift a star's HB location towards higher effective tem- 
peratures and luminosities, depending on the amount 
of rotation. This scenario predicts a sloped HB with 
the shift towards higher luminosities (and hence lower 
gravities) i ncrea sing with effective temperature (Rood 
& Crocker |l989|) . 

Helium-mixing scenario: The observed abundance 
variations in globular cluster RGB stars show that 
these stars are able to mix nuclearly processed mate- 
rial from the vicinity of the hydrogen-burning shell out 
to the surface (e.g. Kraft 1994 ). In particular, the ob- 
served Al enhancements indicate that the mixing can 
penetrate into the H-shell, thus leading to the dredge- 
up of helium. The resulting increase in the envelope 
helium abundance produces an HB morphology that 
slopes upward towards brighter luminosities with in- 
creasing effective temperature (Sweigart 1997 ; SC98). 
The shift towards lower than canonical gravities due to 
helium mixing reaches a maximum between 10,000 K 
and 20,000 K. No shift is predicted for the hottest HB 
stars (>20,000 K), where the H-shell is inactive and 
the luminosity is therefore unaffected by helium mix- 



ing (Sweigart 1997). This scenario might also help to 



explain the low gravities of the blue HB stars found in 
several metal-poor globular clusters (Moehler 1999] and 
references therein), although Grundahl et al. (1999) 



have recently provided some caveats about it. They do 
note, however, that "helium mixing stands out as the 
best candidate to explain the anomalous HB morphol- 
ogy of the metal-rich globular clusters NGC 6388 and 
NGC 6441". 

These scenarios make different predictions for the sur- 
face gravities of HB stars which can be tested observation- 
ally. We started a program to obtain spectroscopic obser- 
vations with this goal in mind and report here on our first 
results. In Sect. |2] we describe our observations and the 
employed data reduction techniques; in Sect. [§] the pro- 
cedure adopted to derive the atmospheric parameters of 
the programme stars is outlined. Finally, we discuss our 
results and their consequences in Sect. [|. 

2. Observations and Data Reduction 

We selected our targets from the HST Wide Field and 
Planetary Camera-2 (WFPC2) images of Rich et al. 
( |1997| ) which were convolved with a "seeing" of 1" in or- 
der to allow the selection of stars suitable for ground-based 
spectroscopic observations. We took spectra of four stars 
in NGC 6388 and three stars in NGC 6441 (see Table 0). 
Unfortunately, due to rather mediocre weather conditions 
we could not observe more stars during our New Technol- 
ogy Telescope (NTT) observing run. 



We used only the blue channel of the ESO Multi-Mode 
Instrument (EMMI) at the NTT, since the use of the 
two-channel mode leads to destructive interference around 
the Hg line. We obtained medium resolution spectra to 
measure Balmer line profiles and helium line equivalent 
widths. We used grating #4 (72 Amm -1 ) and a slit width 
of l'.'O, keeping the slit at parallactic angle for all observa- 
tions. Seeing values varied between 0'.'8 and 1'.'5. The CCD 
was a Tek 1024 x 1024 chip with (24 /mi) 2 pixels, a read- 
out-noise of 5.7 e~ and a gain of 2.84 e~ADU~ 1 (where 
as usual "ADU" stands for analog-to-digital- units). 

For calibration purposes we observed each night ten 
bias frames and ten dome flat-fields with a mean expo- 
sure level of about 10,000 ADU each. In the beginning of 
each night we took a long HeAr exposure for wavelength 
calibration. Due to the long exposure times for the Ar 
spectrum we only obtained He calibration spectra dur- 
ing the night (before and after each science observation) 
and determined the offsets relative to the HeAr spectrum 
(from which the dispersion relation was derived) by corre- 
lating the two spectra. For the second night (and part of 
the third night) we could not obtain any wavelength cali- 
bration spectra due to technical problems. To correct any 
spectral distortions we used HeAr spectra from the first 
night for those data. The FWHM of the HeAr lines was 
measured to be 5.67±0.2 A and was used as instrumental 
resolution. As flux standard stars we used LTT 7987 and 
EG 274. 

We averaged the bias frames of the three nights and 
used their mean value instead of the whole frames as there 
was no spatial or temporal variation detectable. To correct 
the electronic offset we adjusted the mean bias by the 
difference between the mean overscan value of the science 
frame and that of the bias frame. The dark currents were 
determined from several long dark frames and turned out 
to be negligible (5.6 ± 4.2 ^). 

The flat fields were averaged separately for each night, 
since we detected a slight variation in the fringe patterns 
of the flat fields from one night to the next (below 5%). 
We averaged the flat fields and then determined the spec- 
tral energy distribution of the flat field lamp by averaging 
the flat fields along the spatial axis. This one-dimensional 
"flat field spectrum" was then heavily smoothed and used 
afterwards to normalize the dome flats along the disper- 
sion axis. 

For the wavelength calibration we fitted 3 rd -order poly- 
nomials to the dispersion relations of the HeAr spec- 
tra. They were adjusted for eventual offsets by cross- 
correlating them with the He spectra belonging to the 
science frame. We rebinned the frames two-dimensionally 
to constant wavelength steps. Before the sky fit the frames 
were smoothed along the spatial axis to erase cosmics in 
the background. To determine the sky background we had 
to find regions without any stellar spectra, which were 
sometimes not close to the place of the object's spectrum. 
Nevertheless the flat field correction and wavelength cal- 
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Cluster 


Star 


02000 


#2000 


V 


B-V 












[mag] 


[mag] 


NGC 6388 


WF2- 


-42 


17 h 36 m 2r!4 


-44°42'56" 


17.440 


+0.334 


NGC 6388 


WF3- 


-14 


17 h 36 m 23 i !5 


-44°43'11" 


16.914 


+0.333 


NGC 6388 


WF3- 


-15 


17 h 36 m 24 i :i 


-44°43'15" 


16.896 


+0.432 


NGC 6388 


WF4- 


-20 


17 h 36 m 2ffi9 


-44°45'10" 


17.183 


+0.341 


NGC 6441 


WF2- 


-24 


17 h 50 m 08:5 


-37°04'09" 


17.567 


+0.386 


NGC 6441 


WF3- 


-16 


17 h 50 m 06 , !2 


-37°02'50" 


17.683 


+0.558 


NGC 6441 


WF3- 


-17 


17 h 50 m 04!5 


-37°02'20" 


17.643 


+0.499 



ibration turned out to be good enough that a linear fit 
to the spatial distribution of the sky light allowed the 
sky background at the object's position to be reproduced 
with sufficient accuracy. This means in our case that af- 
ter the fitted sky background was subtracted from the un- 
smoothed frame we do not see any absorption lines caused 
by the predominantly red stars of the clusters. The sky- 
subtracted spectra were extracted using Home's ( 1986| ) al- 
gorithm as implemented in MIDAS (Munich Image Data 
Analysis System). 

Finally the spectra were corrected for atmospheric ex- 



tinction using the data of Tug (1977). The data for the 
flux standard stars were taken from Hamuy et al. (1992) 
and the response curves were fitted by splines. 



3. Atmospheric Parameters 

To derive effective temperatures, surface gravities and he- 
lium abundances we fitted the observed Balmer and he- 
lium lines with stellar model atmospheres. Beforehand we 
corrected the observed spectra for wavelength shifts in- 
troduced by radial velocities and the lack of appropriate 
wavelength calibration frames (see Sect. ||), derived from 
the positions of the Balmer lines. The individual spectra 
for each star were then co-added and normalized by eye 
and are plotted in Fig. |l|. For the normalization we as- 
sumed that the observed noise is mostly caused by photon 
noise. However, to account for the possibility of faint un- 
resolved metal lines we tried to place the continuum line 
at the upper region of the scatter. 

To establish the best fi t we used the routi nes de veloped 
by Bergeron et al. ( |l992[) and Saffer et al. (|l994|) , which 



employ a \ test. In addition the fit program normalizes 
model spectra and observed spectra using the same points 
for the continuum definition. The a necessary for the cal- 
culation of x 2 is estimated from the noise in the continuum 
regions of the spectra. 

We computed model atmospheres using ATLAS9 (Ku- 
priv. comm.) and used Michael Lemke's ver- 
of the LINFOR program (developed originally by 
Holwcger, Steffen, and Steenbock at Kiel University) to 
compute a grid of theoretical spectra which include the 



ruc z [1991 
sionP] of t 




4000 



4500 



3000 



a r a ~ 



1 [http: / /www. stcrnwartc.uni-crlangcn.de/~ai26/linfit /linfor.html 



Fig. 1. Normalized spectra for the programme stars 



Balmer lines H a to H22 and He I lines. The grid covered 
the range 7500 K < T cS < 20,000 K, 2.5 < log g < 5.0, 
-3.0 < log^f < -1.0, at metallicities of [M/H] = -0.75 
and [M/H] = -0.50. To check the effects of metallic- 
ity we fitted the spectra at log^ 1 = —1.0 using model 
spectra with [M/H] = -0.50 and [M/H] = -0.75. The 
fit results for the two metallicities were almost identi- 
cal: T e ff and log g, respectively, increased by at most 1% 
and 0.03 dex from [M/H] = -0.50 to -0.75 (except for 
NGC 6441 WF3-16, where AT off and A log g were 1.6% 
and 0.07 dex, respectively). All further analysis was there- 
fore performed for a fixed metallicity of [M/H] = —0.50. 
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As none of the programme stars showed any He ab- 
sorption lines we first fitted all spectra at two fixed helium 
abundances, log = —1.0 and —2.0, using the Balmer 
lines Hp to Hio (except H e to avoid contamination by the 
Call H line). When defining the fit regions we took care 
to extend them as far as possible while avoiding any pre- 
dicted strong He I lines (AA 4026 A, 4388 A, 4921 A) which 
might be present at the noise level although undetected 
by eye. Such lines would distort the fit: Fitting regions 
of strong He I lines while keeping the helium abundance 
fixed will result in erroneous results for T e g and log g if 
the helium abundance of the model atmospheres is not 
identical with that of the star. The results are listed in 
Table || 

As the stars are rather cool it is not obvious from the 
spectra on which side of the Balmer maximum they lie. 
However, the strength of the Call K line can be used 
to distinguish between the "hot" and the "cool" solu- 
tion: the observed equivalent width of this line has stellar 
and interstellar contributions. As the clusters are highly 
reddened - E(B - V) = 0.40 and 0.44 for NGC 6388 
and NGC 6441, respectively (Harris 1996) - the observed 



equivalent width has to be significantly larger than the one 
predicted by the model atmosphere for [M/H] = —0.5. In 
Table || we list the atmospheric parameters for all stars to- 
gether with the measured and predicted equivalent widths 
for the Call K line. The comparison between observed and 
predicted values places all ambiguous stars on the hot side 
of the Balmer maximum. 

For those stars that are hotter than 9,500 K we fi- 
nally fitted the Balmer lines mentioned above together 
with those regions of the spectra where the He I lines 
AA 4026 A, 4388 A, 4471 A, and 4921 A are expected. 
This way we want to put at least an upper limit to the 
He abundance even though no He lines are visible. As can 
be seen from Table | this upper limit tends to be subsolar 
for the hotter stars, in agreement with the results for blue 
HB stars in metal-poor clusters (e.g. Moehler et al. 1997, 
19994 Behr et al. [I999p. 



3.1. Possible errors and discrepancies 

1. Error estimates: The fit program gives la errors 
derived from A X 2 = 2.71 (T eff , log g) resp. 3.53 (T eff , 
log g, log t^). However, since x 2 is not close to 1 these 
errors will most likely underestimate the true errors. 
As the errors causing the large x 2 values are most prob- 
ably systematic (see below) they are hard to quantify. 
We decided to get an estimate of their size by assuming 
that the large x 2 is solely caused by noise. Increasing 
a until x 2 = 1 then yields new formal errors for Ax 2 = 
2.71 and 3.53, respectively, which are given in Table pi 
We therefore adopt a general overall error in log g of 
±0.2 dex and in Tgg of 5%. Fig. || shows theoretical 
line profiles for the parameters derived for NGC 6441 



WF3-17 and model atmosphere spectra for Alog g= 
±0.2 dex. 

2. Large x 2 : As can be seen from Table || the reduced 
X 2 is far from 1. We believe that these bad x 2 va l- 
ues are due to the low resolution of the data (5.7 A): 
At T cff = 11,000 K and log g = 4 the FWHM of the 
theoretical Balmer lines are around 12 A. Thus the in- 
strumental profile makes up a considerable part of the 
observed line profile. Therefore any deviations of the 
instrumental profile from a Gaussian (which is used to 
convolve the model atmospheres) will lead to a bad 
fit. Fortunately we could compare the effects for one 
blue HB star in NGC 6752, for which we have an NTT 
spectrum with the same setup as is used here and a 
spectrum from the ESO 1.52m telescope with areso- 
lution of 2.6 A (for details see Moehler et al. |l999b|) . 
Performing a spectroscopic analysis as described above 
with model atmospheres of JM/H] = —4.5 we find the 
parameters listed in Table ||. While the \ 2 values for 
both fits differ considerably the resulting effective tem- 
peratures and surface gravities are still rather similar. 

3. Continuum placement: Using the spectrum of 
NGC 6441 WF3-47 we checked the effects of a pos- 
sible misplacement of the continuum: We now nor- 
malized the spectrum again, using only the uppermost 
resp. lowermost continuum points. This amounted to 
changes of —1% (upper continuum) resp. ±3% (lower 
continuum) in the normalized continuum level. Tem- 
perature, gravity and x 2 were only slightly affected by 
these changes: In both cases the temperature increased 
by about 50 K, log g by 0.03 dex, and x 2 by about 0.04. 
The effect of a misplaced continuum is mostly erased 
by the renormalization of observed and model spec- 
trum before fitting, using the same continuum points 
for both. 

The effective temperatures from Table || and the B — V 
colours from Table [l] do not correlate well, esp. for the 
stars in NGC 6441. However, as the reddening towards 
these globular clusters is relatively large [E(B — V) « 
0.4 mag] and may vary by a few hundredths of a magnitude 
on small scales (especially in the case of NGC 6444; see 
Sect. 3.1 in Piotto et al . |1997 ; for a more general discussion 
see Heitsch & Richtler 1999), we do not put too much em- 



phasis on this (apparent) inconsistency between spectro- 
scopically derived temperatures and observed (WFPC2) 
colours. 



4. Discussion and Future Work 

The atmospheric parameters derived for our programme 
stars are compared in Fig. || to both the canonical ZAHB 
and the non-canonic al tracks fo r [M/H] = —0.5 from 
Sweigart & Catelan (1998, 199S , helium-mixed resp. ro- 
tation, cf. Sect. fjj). Quite unexpectedly, the studied stars 
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Table 2. Atmospheric parameters and equivalent widths of the Call K line for the programme stars for an assumed 
[M/H] = —0.5. For stars hotter than 9,500 K the first row gives the atmospheric parameters obtained by fitting the 
Balmer lines and the spectral region of the strongest He I lines. For all stars the first two of the four rows at fixed 
helium abundance give the solutions on the hot side of the Balmer maximum for the cited helium abundances, while 
the last two rows give the solutions on the cool side of the Balmer maximum that were rejected because of the strength 
of the Call K line. The errors are la errors adjusted for x 2 > 1 (see also Sect. |||). A colon marks extrapolated values 
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tend to lie preferentially below the canonical ZAHBq This 
behaviour stands in marked contrast to the results for 
blue HB stars in more metal-poor globular clusters, where 



2 The d iffere nce between Fig. ^ in this paper and Fig. 8 in 
Moehler ( 1999 ) results from using fitting ranges for the pre- 
liminary analysis of the line profiles that did not extend far 
enough to include the wings of the Balmer line profiles. 



lower than canonical gravities are normally found (see 



Moehler 1999 and references therein). To verify that the 



high gravities are not due to problems with the low reso- 
lution of the data, we also show results for cool blue HB 
stars in NGC 6752 that were observed du ring th e same run 
with the same resolution (Moehler et al. 1999b ). ft can be 
clearly seen that those stars follow the overall trend seen 
in NGC 6752 to show lower than expected surface grav- 
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Fig. 2. Plots of model atmosphere spectra for T e g = 10,910 K and log g = 4.04 (solid line), 3.84 (upper dashed line), 
4.24 (lower dashed line). The model spectra were convolved by a Gaussian of 5.67 A FWHM to take into account the 
instrumental resolution of the observed data. Given are those ranges of the Balmer lines that were fitted to obtain 
effective temperatures and surface gravities. The asymmetric fit regions are due to the avoidance of He I lines that were 
fitted separately (see Sect. ||). The right panel shows the observed line profiles of NGC 6441 WF3-17 in comparison 



Table 3. Results of spectroscopic analyses for star B 3253 
in NGC 6752 
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ities when fitted with metal-poor model atmospheres. If 
the high gravities found for the stars in NGC 6388 and 
NGC 6441 were caused by systematic errors in the data 
analysis the stars in NGC 6752 should exhibit the same 
effect. 

Considering the estimated errors in T e g and log g the 
positions of NGC 6388 WF2-42, WF3-14, WF4-20, and 
NGC 6441 WF3-17, however, are consistent with the 
canonical ZAHB within 2a. NGC 6388 WF3-15 could be 



a blue HB star already on its way to the asymptotic giant 
branch - as we observed rather bright blue HB stars at the 
top of the blue tail in both clusters we are biased towards 
stars evolving from the ZAHB towards higher luminosities. 
NGC 6441 WF3-16 lies marginally more than 2a below 
the unmixed ZAHB, but has also the lowest S/N of all 
observed stars. NGC 6441 WF2-24 lies 0.7 dex in log g 
below the canonical ZAHB along th e trac k of a 0.195 M Q 
helium white dwarf (Driebe et al. 1998 ). In this case it 
would be a foreground object and most probably mem- 
ber of a binary system having undergone mass transfer, 
as the low mass single star precursors of helium white 
dwarfs evolve too slowly to reach this stage within a Hub- 
ble time. While helium white dwarfs evolve rather fast and 
the probability of finding one is therefore rather low, we 
could not find any other physical explanation for a star 
like NGC 6441 WF2-24. As we are not able to derive ra- 
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Fig. 3. Comparison of the measured gravities and tem- 
peratures for our programme stars in NGC 6388 and 
NGC 6441 against th e HB models for [M/H] = -0.5 from 
Sweigart & Catelan (1999). The solid curve represents a 
canonical zero-age HB (ZAHB) for this metallicity. The 
ZAHB and evolutionary tracks for the sequences which re- 
produce the observed HB slope in these globular clusters 
are indicated by the dashed and dotted curves, respec- 
tively. The upper panel shows the helium-mixed tracks, 
the lower panel shows the rotation tracks. The numbers 
of the stars refer to Table ||. A representative error bar is 
plotted at the location of the star NGC 6441 WF3-17. In 
addition we show cool blue HB stars in NGC 6752, which 
were observed with the same instrumental setup during 
the same run (Moehler et al. 1999b) 



dial velocities (see Sect. ^) we cannot decide whether this 
star is a cluster member or not. 

The derived gravities (except for NGC 6388 WF3- 
15) are significantly larger than those predicted by the 
non-canonical tracks that reproduce the upward sloping 
HB's in the colour-magnitude diagrams of NGC 6388 and 
NGC 6441 (see Sect. |l|). The discrepancy would increase if 
the atmospheric abundances of the programme stars (and 
thereby their line profiles) were affec ted b y the radiative 
levitation of metals (Grundahl et al. 1999). As shown by 



parameters of blue HB stars in NGC 6752 to lower tem- 
peratures and/or higher gravities. However, the moderate 
resolution and S /N of the data discussed here do not allow 
abundances of the heavy elements to be estimated. 

We do not have an explanation for this surprising re- 
sult. All the scenarios laid out by SC98 and discussed in 
Sect. [I] predict anomalously low gravities for cluster blue 
HB stars within the temperature range of our programme 
stars. As far as we are aware, there are no alternative 
models capable of accounting for the sloped HB's seen 
in NGC 6388 and NGC 6441 without producing anoma- 
lously bright HB stars and hence low gravities. In fact, 
recent analyses of RR Lyrae var iables in NGC 638 8 and 
NGC 6441 (SC98; Layden et al. |l999t Pritzl et al. |l999[ ) 



strongly indicate that the HB's of these globular clusters 
are substantially brighter than canonical models would 
predict. Thus we face a conundrum: the non-canonical 
models which explain the upward sloping HB's in these 
globular clusters are inconsistent with the derived gravi- 
ties. 

It is clear that further spectroscopic observations of 
larger (statistically more significant) samples of blue HB 
stars in NGC 6388 and NGC 6441 are urgently needed to 
verify the present results. 
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